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Abstract 

The effect of surface roughness on the Seebeck coefficient in the sub-50-nm scale silicon ultra thin films is 
investigated theoretically using nonequilibrium Green's function formalism. For systematic studies, the surface 
roughness is modelled by varying thickness periodically with square wave profile characterized by two parameters: 
amplitude {Aq) and wavelength (A). Since high Seebeck coefficient is obtained if the temperature difference 
between the ends of device produces higher currents and higher induced voltages, we investigate how the 
generated current and induced voltage is affected with increasing Aq and A. The theoretical investigations show 
that pseudoperiodicity of the device structure gives rise to two effects: firstly the threshold energy at which the 
transmission of current starts is shifted towards higher energy sides and secondly transmission spectra of current 
possess pseudobands and pseudogaps. The width of the pseudobands and their occupancies determine the total 
generated current. It is found that current decreases with increasing Aq but shows a complicated trend with A. The 
trends of threshold energy determine the trends of Seebeck voltage with roughness parameters. The increase in 
threshold energy makes the current flow in higher energy levels. Thus, the Seebeck voltage, i.e. voltage required to 
nullify this current, increases. Increase in Seebeck voltage results in increase in Seebeck coefficient. We find that 
threshold energy increases with increasing Aq and frequency Hence, Seebeck voltage and Seebeck coefficient 
increase vice versa. It is observed that Seebeck coefficient is tuneable with surface roughness parameters. 



1. Introduction 

It is reported that when the two ends of the thermoelec- 
tric material are kept at different temperatures, the vol- 
tage is generated [1]. This effect is called as the Seebeck 
effect. This effect has useful applications as it shows that 
heat can be converted to electricity and vice versa. In the 
last few years due to quest for finding, the highly efficient 
ways to produce and use energy there has been an 
increased interest in the thermoelectric performance of a 
device [2]. Moreover, thermoelectric energy conversion 
paves the way for utilization of waste heat generated in 
the devices into useful electric power. The extraction of 
waste heat is an important issue, especially nowadays, 
when the devices are reaching the sub-50-nm scale. 
Thus, intensive research has been reported for the effect 
of dimensionality on the thermoelectric phenomenon 
[3-5]. As the dimensions reduce down to sub-50-nm 
scale, the surface roughness plays a crucial role in 
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determining the performance of the device. In this scale, 
the spatial extent of variations due to roughness is com- 
parable to the device dimensions [6] . Hence, the investi- 
gation on the effect of surface roughness is important 
and crucial [7-9]. In this article, we have made a systema- 
tic study of the surface roughness effects on Seebeck 
coefficient. The surface roughness has been modelled by 
varying the thickness of the film periodically with a 
square wave profile as shown in Figure 1. The roughness 
is characterized by two parameters: amplitude (Aq) and 
wavelength (A) of the square wave [10]. These parameters 
are analogous to the root mean square roughness and the 
roughness autocorrelation length which are commonly 
used to determine the surface roughness morphology 
[11]. We study the effect of roughness with increasing Aq 
and A on the Seebeck effect. In other words, we investi- 
gate how does the induced voltage, Vseebecio for a fixed 
temperature difference between the two ends of a device, 
is affected due to increasing Aq and A of the surface 
roughness of the device. To clearly bring the role of sur- 
face roughness the electron transport is considered to be 
coherent inside the device. The phase coherent transport 
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Figure 1 The schematic diagram of simulated 2D Si structure 

Length of the channel is L^, average thickness is Tb and roughness 
parameters are: amplitude {Aq) & periodicity (k). The bias is applied 
to nullif/ the current due to temperature gradient. 



is possible in the devices with the dimensions in the sub- 
50-nm scale [12]. 

When the two ends of the device are subjected to a 
temperature difference, initially, there is a net flow of 
current I^et from one end of the device to the other as 
shown in Figure 2a. This is so because, due to tempera- 
ture difference, there is a different distribution of elec- 
trons in different energy states at the two ends. Hence, 
the number of electrons moving from the cold end to 
the hot end in an energy level 's' is not the same as that 
of from hot end to cold end. This results in net flow of 
electrons in an energy level 'e'. The total current I^et is 
obtained from adding the contributions from all the 
energy levels. This net flow of current results in an 
accumulation of electrons at one end which builds an 
electric field opposing the incoming electrons as shown 
in Figure 2b. Further accumulation of charge stops 
when no more electrons are able to make it due to the 
building of electric field. The voltage developed due to 
the electric field is called Seebeck voltage Vseebeck- See- 
beck coefficient is obtained by dividing the Seebeck vol- 
tage by temperature difference [13]. As shown in Figure 
2c, Seebeck voltage can be measured by calculating the 
amount of voltage required 'VappHed' to nullify the cur- 
rent '/net' due to temperature gradient. In this article, we 
focus on two important issues: firstly we investigate how 



does, for the same temperature difference between the 
two ends, /net changes with Aq and A. Secondly, it has 
been observed that for some roughness parameters, 
though the /net decreases but the voltage applied to nul- 
lify it, Vseebeck> increases. We investigate the reasons 
behind this behaviour. 

2. Theory 

The expression for the current has been calculated using 
nonequilibrium Green's function formalism (NEGF) 
[14]. For the devices in the sub-50-nm scale quantum 
effects become important and quantum transport model 
is required for description of transport phenomenon. In 
this regards, NEGF formalism, in which quantum effects 
are inherent, provides an efficient framework to model 
the electron transport in thermoelectric devices. A brief 
description is given below and detailed development can 
be found in [15]. 

The isolated channel and its energy levels are 
described by Schrodinger equation using Hamiltonian 
H, potential U and energy eigenvalues of the electron e^^ 



(1) 



But when the channel gets connected to the contacts 
it becomes an open system and the electrons can come 
and go out of the channel. To represent such an open 
system, two extra terms get introduced in the Schrodin- 
ger equation. The Schrodinger equation becomes 



{EI - H - j:^ - j:^) {if] = {s] 



(2) 



El and Z2, called self-energy, represent the coupling 
between the channel and the two contacts. Z's are a 
complex quantity. Its imaginary part signifies that the 
electrons can decay out of the channel. 



Decay rate 
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The real part of Z's just lifts up the energy eigenvalues 
of the channel. The second term that gets introduced in 
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Figure 2 Electron flow for different configurations of device, (a) Electrons flow from hot end to cold end due to temperature gradient, (b) 
An induced emf is generated due to temperature gradient (c) Voltage is applied to nullify the current due to temperature gradient. The voltage 
at which current becomes zero is called Seebeck voltage. 
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the Schrodinger equation is {S}, it is called the source 
term and is like an actual source driving the system. 
Solution of the Schrodinger equation is given by 

iA = [G] [S] 
where 

[G] = (£/-H-Ei-E2r' 



- EI-H-E- EI-H-Re (E)-nm^(E) ^ ' 

[G] is called the Green's function of the system. In the 
energy domain, Green function gives the energy eigen- 
values for the eigenstates that are occupied in response 
to the applied impulse. 

Density of states at energy's' = A (£) = i (G - G*) = GFG* (5) 
Asr = Ti + 

A[E) = GViG^ + Gr2G+ 

A(£)=Ai(£)+A2 (£) 

Ai(E) and A2{E) represent density of states due to 
coupling of the channel with contacts 1 and 2, respec- 
tively. The probability of filling of these states is deter- 
mined by the Fermi Dirac distribution function, y*(£) 
and^(£), of the contacts 1 and 2, which are at tempera- 
ture T and T+AT respectively. 

Carrier concentration is given by: 

n(£) =Ai (£)/!(£)+ A2(£)f(£) 
= G" (£) 

The net inflow of current in channel at contact 1 is 
the same as the net outflow of current at contact 2. Net 
current flow in channel at contact 1 is the difference 
between the inflow and outflow of current at contact 1. 

/ (£) = Ti (£) /I (£) A (£) - Ti (£) G" (£) 

= Ti (£) G (£) (£) G* (£) {f (£) - f (£)) 
= T(£) (fi (£)-/2(£)) 

If/(£) = 1 and/(£) = 0 then £(£) represents the 
maximum current flow in the channel. For the electron 
transport in ;«-direction with confinement present in z- 
direction, the expression for 1(E), after integrating over 
the energies along j/-direction, becomes 

/ (£^) = T (£,) (fi^ (£^) - (£,)) (7) 



Total current = I = J2fl (Ex)dEx (g) 

where is the subband index. 
3. Results and discussion 

The information about how does Ineo i e- the current 
due to the temperature gradient, and, Vseebecio i e- the 
voltage required to nullify it, depend on Aq and A, is 
obtained by plotting the current voltage {I-V) character- 
istics for devices with different surface roughness para- 
meters. Figure 3 gives the /-V^ characteristics for (a) 
different amplitudes: Aq's = 0.1 nm, 0.3 nm and 0.5 nm 
and A = 2.5 nm and (b) different wavelengths: A's = 2.5 
nm, 3.3 nm and 5 nm and Aq = 0.5 nm. The plot for 
smooth surface is given in Figure 3a for comparison. 
The temperature gradient of A£ = 20 K is applied 
across the device. The graphs are obtained for the 
devices in the sub-50-nm scale by considering the para- 
meters: length of the channel Lc = 10 nm and average 
thickness T/, = 2 nm. As /net is the current generated 
only due to temperature gradient with Vappiij.^ = 0, it is 
given by the intercept with the j-axis. It is observed 
from Figure 3a, b that /net decreases with increasing Aq 
but with A there is no definite trend. Vgeebeck is the vol- 
tage applied to make the current due to temperature 
gradient equal to zero and is obtained by the intercept 
with the X-axis. It is observed from Figure 3a, b that 
Vseebeck increases with increasing Aq and 1/A. 

The above trends of /net and Vseebeck with the rough- 
ness Aq and 1/A can be understood by examining the 
current flowing through each energy level /(£). As seen 
from Equation 7, 1(E) is given as the product of trans- 
mission £(£) and (^/[g''^ (£) - /f^' (£)^ , i.e. difference in 

the occupancy of cold /i™''' (£) and hot junctions 
fiQ (£) • The transmission £(£) gives the maximum cur- 
rent through the energy level '£' and depends on the 
density of states A{E) and velocity, i.e. rate of flow of 
electrons in and out of that energy level. As such the 
geometry of the device does not affect the occupancy 

and hence (/[°''^ (£) - °' (£) j but it affects the trans- 
mission £(£) of electrons from energy level £ through 
density of states A{E). Thus, the trends of I^t and Vsee- 
beck with Aq and A can be understood by examining how 
the transmission gets affected by roughness. 

3.1. Effect of roughness on transmission spectra 

Figure 4 gives the plots of £(£) for (a) roughness wave- 
length: A's = 2.5 nm, 3.3 nm and 5 nm with Aq = 0.5 
nm and (b) roughness amplitude: Aq's = 0.1 nm, 0.3 nm 
and 0.5 nm with A = 2.5 nm. The following trends are 
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Figure 3 Current-voltage characteristics for the 2D Si structure. U= 10 nm, Tt = 2 nm, transport mass rVx = 0.19 mo, quantization mass 
= 0.91 mo and AT = 20 K. For (a) 1 = 2.5 nm Mo's = 0.1 nm, 0.3 nm, 0.5 nm and for (b) Ao = 0.5 nm &;i's = 2.5 nm, 3.3 nm, 5 nm. Also in (a) l-V 
cliaracteristics for smooth surface is sliown for comparison. 



observed: (i) threshold energy, i.e. the 'onset' energy at 
which the transmission starts, increases with increasing 
^0 and 1/A. (ii) There are regions of psuedobands and 
pseudogaps. The width of the pseudobands decreases 
with increasing Aq and A. Within the pseudobands, 
there are peaks in transmission at some energy. The 
phenomenon of increase in threshold energy has also 
been observed experimentally in quantum wells with 
thickness < 4 nm [16]. We would like to mention here 
that it has been observed only for thicknesses below 4 
nm because surface roughness effects become promi- 
nent for ultra thin films. The criterion for the thickness 



below which surface roughness effects become promi- 
nent is discussed in [16]. So, our discussions can be 
experimentally validated only for ultra thin films. We 
will analyse the implications of the above trends on /net 
and Kseebeck but before that we discuss in brief the pre- 
sence of the above trends. 

The increase in threshold energy occurs due to the 
varying thickness seen by the electron while crossing the 
channel. The regions of the channel where the thickness 
is large, the confinement is small and the regions where 
the thickness is small, the confinement is high as seen 
from Figure 1. Since the lowest value of e^, which 




0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 
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Figure 4 Energy resolved transmission spectra at zero applied voltage. For (a) Aq = 0.5 nm &A's = 2 5 nm, 3 3 nm, 5 nm and for (b) A = 

2.5 nm &4o's = 0.1 nm, 0.3 nm, 0.5 nm. Temperature gradient AT = 20 K for both (a, b). 
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Figure 5 Suband energy profile along the length of the device 

Roughness amplitude (/\o) = 0.5 nm and periodicity (A) = 2.5 nm. 
The regions of channel at which confinement is high, subbands 
energies are high and vice versa. 



corresponds to the subband energy, depends on the con- 
finement it alternates with the same periodicity as the 
thickness as shown in Figure 5. The value of e^, i.e. sub- 
band energy is large for high confinement and vice versa. 
The subband energy corresponds to the potential energy 
seen by the electron. The electron crossing the channel 
experiences this varying potential energy landscape 
which was not the case when the surface was smooth. As 
the electron crossing the channel moves from the region 
of low potential energy to the high potential energy it 
feels an additional confinement effect which was not pre- 
sent for smooth surface. This additional confinement 
effect causes the increase in lowest energy which the 
electron can take. Thus, the threshold energy, which cor- 
responds to the lowest energy eigenvalue, increases. The 
confinement effect will be more for increasing Aq as the 
electron while moving from low potential energy to high 
potential energy region sees higher barrier heights. Thus, 
the lowest energy eigenvalue, i.e. threshold energy, 
increases with increasing Aq. The increase in threshold 
energy with increasing 1/A can be explained on the simi- 
lar lines. The longer A with the same Aq corresponds to 
the wider regions of the low potential energy with the 
same barrier height. This means the confinement effect 
reduces with increasing A. Hence, the threshold energy 
decreases with increasing A. In other words threshold 
energy increases with increasing 1/A. 

The presence of pseudobands and pseudogaps, a fea- 
ture of only rough surfaces, is also due to alternate 
regions of low and high potential energy seen by the 
electron while crossing the channel. The energy eigenva- 
lues for an electron in such a periodic potential energy 
profile are given by Kronig Penny model. The solution 
of the Schrodinger equation for such a potential energy 
profile shows the presence of energy bands and energy 
gaps. The investigations on the bandwidth (BW) of the 



pseudobands with respect to roughness parameters is 
important as the width of the pseudoband determines 
the area under the transmission curve and play an 
important role in determining the total current /^et. It is 
seen from Figure 4a, b that BW of the pseudobands 
decreases with increasing Aq and A. These trends of BW 
on roughness parameters can be understood by looking, 
again, at the Kronig Penny model. According to Kronig 
Penny model, if the barrier heights and widths of the 
periodic potential profile are high, then the BW's are 
small. This is so because higher barrier heights and 
widths give smaller tunnelling probability. As the tun- 
nelling probability becomes smaller the BW reduces. 
Along the same lines increasing Aq, which corresponds 
to higher barrier heights, and increasing A, which corre- 
sponds to larger barrier widths of the pseudoperiodic 
potential profile, would give rise to smaller BW's. Thus 
BW decreases with increasing Aq and A. 

3.2. Effect of roughness on total current /net 

As /net represents the current only due to temperature 
gradient with VappUed = 0, it is given by the intercept with 
the y-axis of /-V characteristics. It is observed from 
Figure 3a that /net decreases with increasing amplitude 
but Figure 3b shows that there is no definite trend 
with A. These features can be explained by considering 
that total current depends on the bandwidth of the pseu- 
dobands and their occupancy. For increasing Aq as seen 
from Figure 4b, threshold energy increases and BW of 
the pseudobands decreases. The smaller bandwidths 
along with their presence at higher energy sides result in 
the decrease in occupancy of these bands. This, in turn, 
decreases the total no. of carriers contributing to the cur- 
rent. Thus, /net decreases with increasing Aq. The trends 
of /net with A are complicated. As seen from Figure 4a 
though the pseudobands are present at lower energy 
sides with increasing A but their BW decreases. Thus, on 
one hand their presence at lower energy sides increases 
the occupancy on the other hand smaller bandwidths 
decrease the total occupancy which is obtained by sum- 
ming the occupancy from all energy levels. These two 
competing features complicate the trends of /net with A. 
Figure 4a shows that BW for A = 3.3 nm is smaller than 
A = 2.5 nm but since the band is present at lower energy 
sides its occupancy is more. This results in more current 
for A = 3.3 nm than for A = 2.5 nm. For A = 5 nm the 
BW is so small that the current is the least. 

3.3. Effect of roughness on Seebeck voltage l/seebeck and 
Seebeck coefficient 

As already mentioned Vseebeck increases with increasing Aq 
and 1/A. This trend can be understood by examining the 
circumstances for which I^^tai (= ^net) becomes zero, /total is 
the sum of current from all energy levels, i.e. / = jI{E)dE. 



Kumar ef al. Nanoscale Research Letters 2012, 7:169 
http://www.nanoscalereslett.eom/content/7/1/169 



Page 6 of 8 



0.8 



0.6 



0.4 



0.2 



0 
0.1 



(a) 








- X=2.5 nm 














I 
















J 


















J 









0.2 



0.3 



0.4 



'e 

< 



0.8 
0.6 
0.4 
0.2 
0 



-0.2 
0.1 



0.2 



0.3 



0.4 





0.2 




0.15 




0.1 




0.05 


'e 




< 


0 


rrent 


-0.05 


u 


-0.1 




-0.15 




-0.2 




-0.25 



0.8 



0.6 



ra 0.4 



0.2 



0.5 



(b) 1 
































\ null 











'e 

< 



0.5 





































. ^null 

If 










— 





































'e 

< 



0 
0.1 



(d) 

1 L 


^=3.3 nm 





























V 







0.2 



0.3 



0.4 



0.5 



1 

0.8 
0.6 
0.4 
0.2 
0 

-0.2 



-0.4 
0.1 

0.4 
0.3 
0.2 
0.1 
0 

-0.1 
-0.2 
-0.3 



. (e) 


































^null 























0.2 



0.3 



0.4 



0.5 



0.1 



-0.4 

0.2 0.3 0.4 0.5 0.1 0.2 0.3 0.4 

Energy (eV) Energy (eV) 

Figure 6 Energy resolved transmission, Fermi Dirac occupancy and current for A = 2.5 nm and 3.3 nm with Ao = 0.5 nm Energy 
resolved transmission is shown in (a, d), difference in the Fermi Dirac occupancy of cold and hot ends is shown in (b, e) and current is shown 
in (c, f). For (a-c) A = 2.5 nm and for (d-f) A = 3.3 nm. fnyn shifts towards higher energy sides with decreasing wavelength. 
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Figure 6c gives I{E) versus E for A = 2.5 nm and Aq = 0.5 
nm. It is seen that below E^^^n the current is +ve and above 
£nuii the current is -ve. Hence, the total current obtained 
by the sum of current from all energy levels is small. We 
need to check when does it become zero. As seen from 
Equation (7), /(£) is product of transmission T(E) and 

(/id'' ~ /id') ' i-^- difference in the occupancy of cold 



/I'p' and hot junctions f^l^^- Since T{E) > 0 for all energies 
therefore current changes sign whenever (fi^'^ — fiu^ 

changes sign as seen from Figure 6b. For those energy 
levels in which the number of electrons at the cold side 
are more than the number of electrons at the hot side, the 

net current flows from cold to hot and ^/jq'^ — /Jd') is 
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+ve and /(£) is +ve. These are the energy levels below fnuii- 
It is vice versa for the energy levels above E^uu- At £nuii> 

(/id'* - /id') = 0- This fnuii is a function of KppHed as 
occupancies of energy states at cold side, f^^'^ and hot 
side, depend on Vappiied- £nuii shifts towards higher 
energy sides with increasing VappUed- We apply voltage to 
shift £nuii to such an extent that total current coming from 
energy levels below fnuii balances the total current coming 
from energy levels above £nuii so that the total current 
becomes zero. 

To understand how does £nuii shifts towards higher 
energy states with increasing Kippiiedi we observe that as 
the threshold energy increases the electrons occupy 
higher energy states and the current contributions come 
from higher energy levels. The net current I^ev which 
flows from hot side to cold side is also made of electrons 
flowing in higher energy levels. We need to apply voltage, 
Vseebeck> to nullify this current I^et- Thus, we need to 
increase the potential of the cold side so that electrons at 
the cold side occupy higher energy states to balance this 
current. Thus, as the threshold energy increases we need 
to apply higher Vseebeck to balance the net current. Also, 
since the occupancy of higher energy states at the cold 
side increases with increasing Vseebeck the £nuii> for which 

(fm'^ ~/id') = 0, shifts towards higher energy sides. 

Since the threshold energy increases with increasing Aq 
and 1/A, hence Seebeck voltage Vgeebeck increases with 
the increasing roughness amplitude and frequency. Since 
Seebeck coefficient is obtained by taking the ratio of Vsee^ 
beck with temperature difference AT, i.e. S = Vjeeback/^^", 
therefore the increase of Vseebeck for fixed AT causes the 
increase in Seebeck coefficient with the similar trends. It 
is seen from Tables 1 and 2 that Seebeck coefficient 
increases with increasing and 1/A. The dependence of 
Seebeck coefficient on roughness parameters suggest that 
it can be tailored by choosing the appropriate roughness 
parameters. 

For the devices if is increased at constant A it is 
observed that though /net decreases but still Vseebeck 
increases. This happens because BW decreases and 
threshold energy increases. Smaller BWs at higher 
energy sides reduce occupancy and hence /net- Though 

Table 1 Variation of Seebeck coefficient for different 



roughness amplitudes 



Ao (nm) 


S (mV/K) 


Smooth 


470 


0.1 


475 


0.3 


537 


0.5 


740 



Roughness wavelength (A) = 2.5 nm is fixed. 



Table 2 Variation of Seebeck coefficient for different 



roughness wavelengths 



X (nm) 


S(nV/K) 


Smooth 


470 


5.0 


545 


3.3 


646 


2.5 


740 


Roughness amplitude {Aq) 


= 0.5 nm is fixed 



the current reduces, it is made up of electrons present 
in higher energy levels and as already explained higher 
voltages are required to nullify the current from higher 
energy states. 

The discussions show that it is the increase in thresh- 
old energy which causes the current to flow in the 
higher energy levels and hence results in increase in 
Seebeck coefficient. This result implies that any physical 
geometry whether periodic or aperiodic which results in 
an increase in threshold energy will show an increase in 
Seebeck coefficient. It also implies that if the threshold 
energy dependence on roughness parameters become 
feeble then Seebeck coefficient will saturate, i.e. it will 
not change much with the change in roughness para- 
meters. The detailed investigations on the above two 
implications are in progress. 

4. Conclusions 

The surface roughness causes the net current Inet> pro- 
duced due to the temperature gradient, to decrease with 
increasing roughness Aq but trends with A are compli- 
cated. The voltage Vseebeck- required to nullify this /net 
shows a definite trend. It increases with increasing Aq 
and 1/A. Increase in Vseebeck results in increase in See- 
beck coefficient. The reasons for the above trends have 
been attributed to the varying periodic thickness seen by 
the electron crossing the channel. The periodically vary- 
ing thickness corresponds to varying confinement and 
hence periodically varying subband energy. The subband 
energy corresponds to the potential energy seen by the 
electron. For the electron moving in a pseudoperiodic 
potential, the transmission spectra shows the following 
features: (1) The threshold energy at which the trans- 
mission of current starts is shifted towards higher 
energy sides. It increases with increasing Aq and 1/A. (2) 
There are regions of pseudobands and pseudogaps. The 
BW of the pseudobands decreases with increasing Aq 
and A. 

The above features of transmission spectra result in 
following trends for /net and Vseebeck^ (1) It is observed 
that BW decreases and shift towards the higher energy 
sides for increasing Aq. Both these features result in 
reduced occupancy of the band and hence total carrier 
concentration. Since /net depends on the total carrier 
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concentration thus Imt reduces with increasing Aq. (2) It 
is observed that pseudobands are present at lower 
energy sides with increasing A but their BW decreases. 
Thus, on one hand the presence of bands on the lower 
energy side increases the occupancy and on the other 
hand the smaller bandwidth reduces the total occu- 
pancy. These two competing features give rise to com- 
plicated trends of /net with A. (3) The trends of 
threshold energy determine the trends of Vseebeck with 
roughness parameters. The increase in threshold energy 
results in the electrons to occupy higher energy states. 
Thus the net current I^et flowing from hot to cold end 
is made up from higher energy levels. To nullify this 
current, the occupancy of higher energy states at the 
cold end is increased by applying higher voltages. Thus 
higher threshold energy requires higher applied voltages, 
Vseebeck> to nullify the current. Since threshold energy 
increases with increasing Aq and hence Vgeebeck and 
Seebeck coefficient increases vice versa. (4) The depen- 
dence of Seebeck coefficient on roughness parameters 
suggest that it can be tailored by choosing the appropri- 
ate roughness parameters. 
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